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A high order series expansion is employed to study the thermodynamical properties of a S = 1/2 
chain coupled to dispersionless phonons. The results are obtained without truncating the phonon 
subspace since the series expansion is performed formally in the overall exchange coupling J. The 
results are used to investigate various parameter regimes, e.g. the adiabatic and antiadiabatic limit 
as well as the intermediate regime which is difficult to investigate by other methods. We find that 
dynamic phonon effects become manifest when more than one thermodynamic quantity is analyzed. 
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I. INTRODUCTION 

In solid state physics, all electronic degrees of freedom 
like charge or spin are coupled to vibrations. Mostly, 
however, such a coupling does not influence the system's 
properties in a decisive way. This is different if the elec- 
tronic degrees of freedom are essentially one dimensional. 
Then the phenomenon of a Peierls transition occurs: the 
system breaks translational invariance spontaneously by 
forming dimers 0, 0, 0> HI El • The interest in a model 
of quantum phonons coupled to one-dimensional spin 
degrees of freedom in particular has been rekindled by 
the discovery of the first inorganic spin-Peierls substance 
CuGeC>3 ■ Single crystals of high quality made investi- 
gations possible that were not possible for the long known 
organic spin-Peierls substances 0. 

Besides the spin-Peierls phenomenon a very strong cou- 
pling of spins and phonons can influence the quantitative 
physics of Mott insulators significantly if the superex- 
changc coupling is small due to geometrical reasons. Ex- 
amples are a 90° angle in the exchange path or a compli- 
cated superexchange process via large ligand groups. In 
these cases, a small change of the geometry implies a cer- 
tain change of the coupling which is very large relative to 
the unchanged coupling. Hence the influence of phonons 
is much larger than usual. Examples for this mecha- 
nism besides CuGeOJEH are SrCu 2 (B0 3 ) 2 [IS El or 
(VO) 2 P 2 7 0ES HIT 

For the above reasons, it is of significant interest to 
provide reliable theoretical predictions for spin systems 
coupled to phonons. It is also clear that acoustic phonons 
will not have a significant influence because they alter the 
exchange pathes relatively weakly. Thus we focus on op- 
tical phonons which have a strong impact on the local ge- 
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ometry so that they can influence the exchange coupling 
significantly. For simplicity, we will consider dispersion- 
less Einstein phonons. The aim of the present work is to 
compute the two fundamental thermodynamic quantities 
of spin systems, namely the magnetic susceptibility and 
the specific heat, and to compare the results in presence 
of a spin-phonon coupling to the results of static spin 
models. This comparison serves as a guideline to experi- 
mental analyses which attempt to identify the signatures 
of spin-phonon couplings. 

The generic spin-phonon model introduced in the fol- 
lowing section cannot be solved analytically. To the au- 
thors' knowledge, there are no analytical exact methods 
to treat extended systems of coupled spins and phonons 
if all energy scales shall be considered. Many approxi- 
mate methods have been applied such as density-matrix 
renormalization |15| , continuous unitary transformations 
EH E3j, exac t diagonalization . linked cluster ex- 
pansion |2dj| . renormalizatio n group |2 l| and quantum 
Monte Carlo (QMC) [3 IH iH 0, l^lig . 

Two limits can be analyzed in more detail. In the 
adiabatic limit u> <C J the spin system is assumed to 
be 'fast' compared to the 'slow' phonon system. Using 
approaches analogous to the ones applied by Pytte p| 
and to the more detailed one by Cross and Fisher Q the 
model in Eq. can be mapped to a statically dimerized 
model. 

The antiadiabatic limit u> 3> J can be handled by an 
appropriate mapping of the starting Hamiltonian to a 
frustrated spin model. Thereby, interactions of larger 
range are induced and the phonon frequency is renor- 
malized E! El El El ES S3 Above a critical frus- 
tration, e.g. a certain next-nearest neighbor interaction, 
the system becomes gapped. This regime is reached for 
large values of the spin-phonon coupling. For small values 
of the spin-phonon coupling the system remains gapless. 
Note that we are dealing with a purely one-dimensional 
problem so that quantum fluctuations may prevent spon- 
taneous symmetry breaking. 

Concerning thermodynamic properties it could be 
shown that the magnetic susceptibility can be fitted by 



2 



a frustrated spin model with temperature independent 
couplings. But this approach fails for increasing values 

Investigations of the regime between the adiabatic and 
the antiadiabatic limit with u> m J are difficult. So far, a 
renormalization group analysis |2ll | and quantum Monte- 
Carlo simulations [13, 123 have been done. 

In the present paper, a spin-phonon model is studied 
using a linked-cluster expansion to derive thermodynami- 
cal properties like the specific heat C and the susceptibil- 
ity x- Zero temperature properties were determined pre- 
viously |2jj. The paper is organized as follows. Sect. UTI 
presents the model, some of its known properties and the 
basic elements of the method employed. In Sect. IIIl| thc 
results are analyzed. They are summarized in Sect. IIVI 
where also open issues are identified. 



II. 



MODEL AND METHOD 



The isotropic spin- 1/2 Heisenberg chain is extended by 
the coupling to local, dispersionless phononic degrees of 
freedom. The Hamilton operator reads 

H = jJ2(l+9(bl + h))SA +1 +LoJ2b\k (la) 

i i 

= H SB +H B . (lb) 

The magnetic exchange coupling is denoted by J, the 
coupling between the phononic subsystem and the mag- 
netic subsystem is given by g J, and the energy of the 
phonons is to. The abbreviations H B and i?sB are used 
in the following. The Hamiltonian in Eq. (JJJ represents 
the so-called bond-coupling model depicted schematically 
in Fig. QJ. The phonons can be viewed to sit between 
the spin sites. They influence only one bond. Such a cou- 
pling can be motivated microscopically [28| for CuGc0 3 . 
But the main reason to choose the coupling as in Eq. 
for our study is its simplicity. 




FIG. 1: Schematic picture of the spin-phonon model as de- 
fined by Eq. QJ. 

In the antiadiabatic limit J/ui — > 0, the critical spin- 
phonon coupling g c for a phase transition from a gap- 
less to a gapped phase is given by g c /uj w 0.4682 using 
the flow equation approach [IS]. Assuming a dimerized 
phase, the model was investigated at T = 20] by a 
linked-cluster expansion which avoided any truncation of 
the phononic Hilbert space. But the starting point was 
the symmetry-broken dimerized phase at zero tempera- 
ture. In the present article, we emphasize the thermody- 
namical aspects of the model without broken symmetry. 



A series expansion about the limit of vanishing J/T is 
performed. The phononic subspace is treated exactly. 
No cutoff in the the phonon subspace is necessary. The 
resulting quantities are given as truncated series in J/T 
with full dependence on the remaining parameters such 
as uj/T. 

First we calculate the partition function Z of the spin- 
phonon system. Then, quantities like the free energy, the 
specific heat or the susceptibility can be derived easily. 
A traditional high temperature series expansion is not 
possible because the expansion in the inverse tempera- 
ture would lead to divergences in the phononic degrees 
of freedom. The limit of infinite temperature is not a 
well-defined starting point for bosons since the phonon 
occupation number diverges. For this reason, we choose 
to perform the formal expansion in the exchange coupling 
J. In this way, the phonon subspace is treated exactly for 
each given temperature. No truncation is necessary and 
the full phonon dynamics is taken into account. To our 
knowledge, this is the first approach of a cluster expan- 
sion about the limit J = at finite temperatures which 
avoids approximations in the phonon subspace. 

We stress that the formal series in J coincides with the 
series for the magnetic subsystem in the inverse tempera- 
ture if we set the spin-phonon coupling to zero: g = 0. So 
it is not surprising that the obtained series bears many 
similarities to a high-temperature series. It is most reli- 
able at high temperatures. The limit to vanishing tem- 
perature is difficult to describe. 

The Hamilton operator from Eq. Q is split into its 
diagonal part Hq = H B and a perturbation V with 

H = H + JV = H B + H SB (2a) 

= c^6i6 4 +J^(l+.9(6l + 6J)S 4 S m (2b) 



Ho V 

The diagonal part Ho is trivially solvable; it describes free 
dispersionless (Einstein) phonons. The perturbation V 
includes the isolated magnetic part and the spin-phonon 
interaction. The standard way to treat such a problem is 
to change to the interaction representation where the off- 
diagonal perturbation governs the non-trivial part of the 
dynamics of the system. In this framework the partition 
function is given as an infinite series in the expansion 
parameter J by 

(3a) 



Z=tr{e-P H } 







= Zo 1 



OO 

„l±-r-^(-J)"/dri 



dr n (V(T 1 )---V(T n )))(3b) 







where the following abbreviations are used. The unper- 
turbed part Ho in (J2J leads to the contribution Zq of the 
partition function 



tr{e-^°} 



n 



(4) 
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with zo = 1/(1 — e /3w ) and the phonon occupation num- 
ber n.j = b\b i . The system size is denoted by N. The 
perturbation V given in the interaction representation as 
V reads 

V( T ) = e THo Ve- rHo (5a) 
= J2 S S :S J+1 (l + ge TH <> (b\ + ft.) e"^°)(5b) 

i 

= Y^SSi+^l + g(bte^ + b i e-^ )) . (5c) 

i 

The angular brackets in Eq. J3J) are an abbreviated no- 
tation for 

(V(n) ■ ■ ■ V(r n )) = J-tr {e-^y(n) • • ■ V(r n )} . (6) 

As can be seen from the above equations the calculations 
for the partition function Z of the magnetic system and 
of the phononic system factorize. In each order of expan- 
sion in J the contribution from the spin system can be 
evaluated separately from the phononic contributions. 

Calculating the partition function in Eq. requires 
repeated integrations over functions of the type 

I(k,l;x n ) = x*e lXn with k G No, ( G Z, i„6l. (7) 

The resulting integrals can be solved exactly with 

x n -i 1 

1^0: J dc n x*e lx " = k\ (-^j 
o 

+±l-» t WJTW&'>"~ (8a) 

1 = 0: J dc n x k n = . (8b) 

o 

These equations allow an iterative evaluation of the mul- 
tiple integrals entering the partition function Z . 

A useful check of the calculations is the limit g = 0. 
This special case yields 

— ^isol. phonons -^isol. spins — -^isol. spins • (9) 

Due to the one-dimensionality of the system under 
study a simple cluster algorithm will be used (for an 
instructive review see Ref. 29]). Therein not only the 
connected clusters are calculated but also the discon- 
nected clusters. The unnecessary calculations of the dis- 
connected cluster are not very costly and we can save 
the book-keeping overhead which would be required oth- 
erwise. The problem of subtracting subclusters occurring 
in the linked cluster expansion algorithm is replaced by 
the evaluation of the lattice constants for a given cluster. 



III. RESULTS 

Here the results for the specific heat C and for the 
susceptibility \ are presented in separate subsections, 
Sects. IhTAI and UTTbI A third subsection Sect. lffl"Cl is 
dedicated to the comparison of the results from static 
spin models and those from the spin-phonon model. 

The bare truncated series provides a first impression of 
the behavior of the considered quantities for various sets 
of parameters. But for quantitative predictions the trun- 
cated series are not sufficient as will be seen in the follow- 
ing. Extrapolation techniques are necessary to improve 
the representation of the results for larger values of J. 
It turns out that the description of higher temperatures 
is easily possible whereas the extrapolation becomes am- 
biguous at low temperatures. We attribute this behavior 
to the fact that the phononic and the magnetic subsystem 
behave at higher temperatures more and more indepen- 
dently. Then our series is essentially a high-temperature 
expansion for the magnetic system which is known to 
work well for higher temperatures. For low temperatures, 
however, possible long-range effects set it which elude our 
approach. 

We benchmark our results relative to QMC data for 
selected sets of parameters [25). 

The results for the spin-phonon model are compared to 
those of pure spin systems. As a reference the exact result 
of the isotropic Heisenberg model [33, OH], [32 is depicted 
in the figures for the susceptibility. The specific heat is 
compared to the specific heat of free phonons plus the 
exactly known result for C(T) for the Heisenberg model 

M 

The coefficients of the series expansion results are avail- 
able upon request. 



A. Specific Heat 

A detailed study of the magnetic properties of the sys- 
tem under consideration also includes the investigation of 
the specific heat. Besides the magnetic susceptibility the 
specific heat is an observable which can be experimen- 
tally easily measured and theoretically easily calculated. 
A direct comparison between theory and experiment is 
often hindered by the fact that the phononic degrees of 
freedom dominate the specific heat. Our model takes 
the influence of a strongly coupled optical phonon into ac- 
count. The additional contribution of acoustic phonons 
is beyond the scope of the present investigation. We as- 
sume that the contribution of the acoustic phonons is 
indeed additive so that it can be accounted for if the 
lattice vibrations are known, e.g. from a dynamic lattice 
model. 

The specific heat C(T) is obtained from the results 
for the free energy per site / obtained from the parti- 
tion function Z using standard relations from statistical 
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physics 



/ = F/N 



1 1 



-InZ 



(10) 



The free energy could be computed to order 13 in J. To 
this end, the contribution of 214 connected and of 470 
disconnected clusters had to be evaluated. To illustrate 
the result the first orders of the free energy series are 
given 

-0f= Iwz) = lnz + J 2 (-P 2 + — — 
HJ N V32 16 w 



3 g 2 P 2 



32 oj 



+J 4 



256 V A' 



2 A3 



-6^ 



(24z 2 - 24z + 6) ^ + (-48z + 24) ^ /3 2 

4 2 \ 

(12 - 24z ) ^ + 48^ ) 0) + 0{ J 5 ) . (11) 



To obtain a detailed insight in the behavior of the spe- 
cific heat Dlog-Pade extrapolations are used. The trun- 
cated series alone are not appropriate for temperatures 
below T < J (not shown) . Fig. [3 depicts the extrapola- 
tion results of the specific heat compared to a superposi- 
tion of the free phonon part of the specific heat given by 



Cf 



(1 - e-^y 



(12) 



and the exactly known result Cs for the isotropic Heisen- 
berg model. Various parameter sets are shown. Simple 
Dlog-Pade extrapolations in J are used for each temper- 
ature point. To this end, the exactly known result of the 
free phonons is subtracted from the series obtained. The 
resulting expression is extrapolated; it starts in second 
order in J. Starting from results for the partition func- 
tion in order 13 in J, the maximum order of extrapolation 
of the remaining specific heat is 10. Ordinary Pade ex- 
trapolations would allow a maximum order of 1 1 , but due 
to the differentiation for the Dlog-Pade extrapolations 
one additional order is lost. After the extrapolation, the 
free phonon contribution is added again. This procedure 
is chosen to deal with series which behave qualitatively 
like the high-temperature series for pure spin systems. 

In Fig. |3 the results for three different sets of param- 
eters to — 1J, g — 0.7 (upper panels), to — 1J, g = 1.5 
(middle panels) and u> = 5J, g = 1 (lower panels) are 
shown. The left plots depict the [n, 2] and the right plots 
the [n,4] extrapolations. In the [n, m] scheme the loga- 
rithmic derivative is approximated by a rational function 
where the polynomial in the numerator has the degree n 
while the one in the denominator has the degree m. 

The [n, 2] and the [n, 4] scheme converge very well. The 
[n, 4] extrapolations are more stable in the low tempera- 
ture regime for both parameter sets. Hence, this scheme 
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FIG. 2: Dlog-Pade extrapolations of the specific heat: var- 
ious extrapolations are compared for three different sets of 
parameters. 



is used in the following. The range of validity can be 
specified by T > 0.15 J as long as the spin-phonon cou- 
pling is smaller or of the same order as lu. For values 
of gJ/u > 1 the extrapolations suffer often from spuri- 
ous poles. In the middle panels the defective extrapo- 
lations are visible. The schemes [5,2], [6,2], and [3,4] 
yield extrapolations which differ significantly from the 
other schemes considered. This is due to spurious poles 
in the integration interval with respect to J. For [6, 2] 
even temperatures above T w J are not described reli- 
ably. 

For small spin-phonon coupling g the specific heat is 
well described by the sum of the phononic Cb and the 
magnetic specific heats Cs of the isolated subsystems, cf. 
first row of panels in Fig. |3 Increasing the spin-phonon 
coupling shifts the specific heat to larger values, cf. sec- 
ond row of panels in Fig. |21 In this temperature regime 
(T > J) already the truncated series yield trustworthy 
results. 

Fixing the spin-phonon coupling g and increasing the 
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phonon frequency lu weakens the visible effect of the 
phonon dynamics, cf. third row of panels in Fig. [21 This 
is to be expected since the phononic subsystems becomes 
more rigid, so that it will not be influenced significantly 
by the spin systems and vice- versa. 



g = 0.66, co = 0.66J 



g = 0.44, co = 0.88J 
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FIG. 3: Dlog-Pade extrapolations of the specific heat com- 
pared to QMC data: the left plot shows data for the values 
g = 0.66, lu = 0.66 J, and the right plot depicts the parameters 
g = 0.44, lu = 0.88J. 



Finally, we compare the extrapolated results from the 
cluster expansion in J to data obtained by QMC. Fig.|3 
depicts two sets of parameters of the [6, 4] scheme and the 
corresponding QMC data. The left panel shows the result 
for g = 0.66, lu = 0.66J. This parameter set implies that 
the system is gapped [l8j . 

In the right panel results are depicted for g = 0.44, 
lu = 0.88J which implies a gapless phase. The statis- 
tical error of the QMC data can be estimated from the 
spread of the data points. The extrapolated series yield 
smooth, continuous results which agree very well with 
the QMC data. For temperatures below T/J « 0.15 the 
dotted lines refer to defective extrapolations which are 
depicted for illustration. These extrapolations are known 
to yield unreliable results beforehand because the extrap- 
olants display spurious poles in the integration interval 
of J. Note that each temperature requires an extrapo- 
lation. The extrapolations at higher temperatures turn 
out to be very stable and not defective so that reliable 
results can be obtained for temperatures T > 0.15 J. For 
lower temperatures, the extrapolations are contaminated 
by spurious poles and should not be trusted. 



B. Susceptibility 

The magnetic susceptibility \ is of special interest since 
it is most easily accessible experimentally. Very often 
static models like the dimerized and/or frustrated spin 
chain can already yield a good description of the suscepti- 
bility of one-dimensional systems. Two cases are conceiv- 
able: cither the appropriate model is indeed a static spin 
model or the static spin model should be seen as effective 



model which incorporates the effects of the spin-phonon 
coupling. Of course, it makes only sense to distinguish 
both cases if there are other experimental probes to dis- 
criminate between them. This will be elucidated in Sect. 

linn 

In the antiadiabatic limit lu > J, detailed investiga- 
tions of the susceptibility were done. The spin-phonon 
chain could be mapped to a frustrated spin chain with 
temperature dependent exchange couplings [l^- The 
corresponding susceptibility could be obtained from a 
high temperature series expansion 33] . It was shown that 
X is only little affected by the temperature dependence of 
the coupling constants. Thus it can be neglected and a 
static model is indeed well justified. This finding agrees 
with previous results |14|. It explains, for instance, why 
the magnetic susceptibility of CuGeC>3 can be fitted so 
surp risingly well by a static frustrated Heisenberg model 

We expect that the mapping to a static spin model 
works less well in the crossover regime to the adiabatic 
limit lu < J. There the effects of the phonon dynamics 
should be more clearly visible. We will return to this 
question in Sect. 1111 CI 

The series expansion of the susceptibility is obtained 
from the previous considerations by incorporating an ex- 
ternal magnetic field. In a first step, the coupling of this 
field to the spins is added to the unperturbed Hamilton 
operator Hq which leads to a modified free energy series 
expansion. In a second step, the susceptibility can be 
derived from the free energy series. 

The unperturbed part Ho of the Hamilton operator 
is extended by a magnetic field term leading to 



H B - hM 



(13) 



with the magnetic field h given in units of g\iB ■ The ad- 
ditional term proportional to the magnetization M com- 
mutes with the free phonon part Hb and with the per- 
turbation V as given in Eq. iffll. Thus, the expression for 
V{t) in Eq. (JSJ is unchanged compared to the calculation 
of the specific heat. The magnetic field is included in Hq. 

The formal expression Eq. @ for the partition function 
is unchanged. From Eq. p3|) we deduce the zeroth order 
contribution Zq which now reads 



tr (e 



tl ( e -f3H Be 0hM} 



( 2cosh | 



flh 
2 



JY 



(14a) 
(14b) 



Taking the logarithm of the partition function Z yields 
^InZ = lnz + In ^2cosh 0^)) ( 15a ) 
+ E (- ^ • • • / *»(V(ti) • • • V(r n )) | (15b) 
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with zq as given in Eq. The angular brackets (■ ■ • ) 
denote the coefficients proportional to N in the trace, 
see Eq. To derive the susceptibility the above equa- 
tion has to be differentiated twice with respect to the 
magnetic field h. Finally h is set to zero 



i d 2 



N 



-\nZ 



(16) 



h=Q 



The susceptibility could be computed up to order 12 
in J. The contribution of 2242 connected and of 2810 
disconnected clusters had to be evaluated. To illustrate 
the result the first orders of the susceptibility series are 
given 
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FIG. 4: Truncated series of the susceptibility: various sets 
of parameters are shown. The exactly known result for the 
Heisenberg model serves as a reference. The left panels illus- 
trate the adiabatic limit whereas the right panels show the 
results in the antiadiabatic limit. 



In Fig. the truncated susceptibility series are de- 
picted for various parameter sets. The energy scales are 
given in units of the magnetic exchange coupling J. The 
exact result of the Heisenberg model serves as a refer- 
ence to illustrate the effects of the additional coupling 
to the phonons. The general feature that the results di- 
verge for temperatures below T < 1.5 J is expected for 
the truncated series. But the qualitative behavior of the 
susceptibility is already discernible. 



The left panels depict the adiabatic regime and the 
right panels illustrate the antiadiabatic limit. The fol- 
lowing conclusions can be drawn from the truncated se- 
ries in the temperature regime T > 1.5 J. Fixing the 
phonon frequency to the overall height of the suscepti- 
bility is lowered for increasing spin-phonon coupling g. 
Such a behavior can be understood in a mean-field treat- 
ment of the spin-phonon coupling. For increasing g the 
effective coupling J c g = J(l + g(tf + b )) increases. This 
shifts the whole susceptibility to lower temperatures com- 
pared to the result of the Heisenberg model with the bare 
magnetic exchange J, cf. discussion in Ref. [l8j. 

For fixed spin-phonon coupling and increasing phonon 
frequency u) this effect becomes less pronounced. For in- 
creasing phonon frequency the magnetic and phononic 
degrees of freedom decouple more and more because the 
phonon system becomes increasingly rigid so that it is 
influenced less and less by the magnetic subsystem. Con- 
comitantly, the spin system is less influenced by the 
phononic subsystem. Thus, for fixed g and uj — > oo the 
magnetic properties are dominated by the antifcrromag- 
netic Heisenberg model. 

The truncated series fail to reproduce the very signifi- 
cant maximum of x(T), cf. Fig.|31 This maximum serves 
as a landmark for many experimental analyses. Thus 
the height and the position of the maximum are of great 
interest. Extrapolations are necessary to extend the rep- 
resentation of the results beyond the radius of conver- 
gence of the truncated series so that the maximum can 
be captured reliably. For the susceptibility the same ex- 
trapolation techniques are applied as for the pure spin 
models described in Ref. j33|. Basically, the truncated 
series is extrapolated using Dlog-Pade approximants in 
an Euler-transformed variable. In contrast to the pure 
spin models the extrapolations of the results of the spin- 
phonon model is not performed in the inverse tempera- 
ture (3, but in the magnetic exchange coupling J. For 
each temperature a separate extrapolation in J has to 
be done. Using standard routines from computer alge- 
bra programs this does not pose more problems than the 
previous extrapolations in (3. 

A more serious restriction concerns the behavior at 
large values of J. For the pure spin systems it was very 
efficient to bias the extrapolations in the inverse temper- 
ature such that the known low temperature behavior was 
captured, see e.g. Ref. |33| . For the spin-phonon system, 
however, much less is known about the excitations at low- 
energies. From the phase diagram depicted in Ref. 0] 
it can be deduced whether the system is gapped or gap- 
less. But the precise value of the gap, let alone the form 
of the dispersion, are not available in the limit J — > oo. 
Note that for the present expansion it is this limit that 
we need to understand, not the limit T — > as in the 
high-temperature series. Hence, we do not attempt to 
bias our expansions in their behavior at large values of 
J. But an improved understanding of the limit J — > oo 
will certainly help to obtain even better extrapolations. 

Fig. [S] shows an overview over the susceptibilities ob- 
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FIG. 5: Dlog-Pade extrapolations of the susceptibility: var- 
ious extrapolation schemes are compared for three different 
sets of parameters. The parameters in the first two rows of 
panels correspond to the gapped regime whereas the lowest 
row depicts results in the gapless regime. 



tained from unbiased Dlog-Pade extrapolations for three 
different sets of parameters. The upper two rows cor- 
respond to results in the gapped regime whereas the pa- 
rameters in the last row correspond to the gapless regime, 
see the phase diagram in Ref . [18| . The left panels depict 
the extrapolations of the form [n, 2] and the right panels 
the extrapolations of the form [n, 4]. 

The position and the height of the maximum are de- 
scribed reliably by both extrapolation schemes [n, 2] and 
[n, 4]. This conclusion is based on the agreement of the 
results for different orders n in Fig. [5j We observe that 
the [n, 4] extrapolations converge better than the [n, 2] 
extrapolations for increasing order. Higher values of m 
in the general [n, m] scheme or odd values of m are likely 
to imply spurious poles. Thus, the [n, 4] scheme is used 
in the following to represent the susceptibility for all pa- 
rameter sets shown in this paper. 

For temperatures T < 0.2 J no results are depicted due 
to spurious poles in the extrapolations in J. The range 



of validity of the [n, 4] extrapolations can be estimated 
to be at least T > 0.25 J independent of the parameter 
sets that we analyzed in our investigations. 
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FIG. 6: Dlog-Pade extrapolations of the susceptibility com- 
pared to QMC data: the left plot shows data for the values 
g = 0.66, uj — 0.66 J (gapped phase), and the right plot de- 
picts data for the values g — 0.44, ui = 0.88 J (gapless phase). 



Finally, we illustrate the reliablity of our extrapola- 
tions by a comparison to QMC data. For two different 
sets of parameters Fig.[S]displays the [7, 4] extrapolations 
of the series and the corresponding QMC data. The val- 
ues g = 0.66, u = 0.66 J in the left panel imply that the 
system is gapped. The values g = 0.44, id = 0.88J in 
the right panel correspond to the gapless regime. The 
agreement between the extrapolated series results and 
the QMC data is very good. For temperatures above 
T/J > 0.2 the results coincide. Below T < 0.2J the 
QMC data and the series data deviate from each other. 
At present, we cannot decide whether the deviations for 
0.1J < T < 0.2 J are due to problems in the QMC simula- 
tion like statistical errors or finite-size effects or whether 
they are due to problems in the series extrapolations. 
Below T « 0.1J, spurious poles occur in the extrapo- 
lated integrands leading to defective extrapolations. For 
illustration, these defective extrapolations are shown as 
dotted lines in Fig. 

We emphasize that the convincing agreement of the 
QMC and the series results in the regime T > 0.2 J sup- 
ports the reliability of these results. In particular, the 
position and the height of the important maximum of 
the susceptibility is described quantitatively. 



C. Comparison to Static Spin Models 

In this paragraph the importance of the spin-phonon 
dynamics shall be highlighted. As mentioned before it is 
often possible to describe the properties of a spin-phonon 
model by a static spin model alone, in particular in the 



8 



antiadiabatic regime. It is to be expected that the ef- 
fects of the dynamic nature of the spin-phonon coupling 
are most prominent in the regime where all energies are of 
similar magnitude. To provide evidence for this hypoth- 
esis we perform the following "theoretical" experiment. 

We start from the extrapolated series data for g = 1 
and ujq = J, see symbols in Fig. El These data shall serve 
as "experimental" input which we will then analyze using 
static spin models. The static spin model considered here 
is the frustrated and dimerized spin chain given by 



(18) 



In this way, we imitate the standard procedure one would 
apply to experimental data. The objective is to see to 
which extent such an analysis yields agreement. In par- 
ticular, we are interested to see where such a description 
remains unsatisfactory. Such an unsatisfactory descrip- 
tion based on static spin models is the signature of the 
dynamic nature of the spin-phonon coupling. 
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FIG. 7: "Theoretical" experiment: The data for g — 1 and 
u>o = J (symbols) is taken as (mock) experimental input. 
It is fitted by the susceptibility of static spin model (solid 
lines). The temperature T max denotes at which temperature 
the susceptibility has its maximum; this value sets the natural 
energy scale of the problem. 

In Fig. [7] the "experimental" susceptibility is ana- 
lyzed by the susceptiblity of static spin models, i.e. the 
isotropic Heisenberg model (HM), the purely frustrated 
spin chain and the dimerized and frustrated chain. Obvi- 
ously, the HM is not appropriate to describe the data de- 
picted by the symbols. But the other two parameter sets 
{a = 0.1, 8 = 0) and (a = 0.12, 8 = 0.05) describe the ex- 
perimental data very well for not too low temperatures, 
i.e. T > 0.5T max . This temperature regime corresponds 
to the the range of temperatures where the extrapolated 
high-temperature series for the dimerized and frustrated 
spin chain are reliable [33|. 

Of course, it is possible to distinguish between different 
models if reliable (experimental) data down to low tem- 
peratures are available. In practice, however, this is often 



not the case since the data at low temperatures can be 
contaminated by impurity effects or other imperfections 
like inclusions of other phases or simply the presence of 
other structural elements in the sample. In such a situ- 
ation, the determination of the appropriate microscopic 
model is difficult and ambiguities are hard to avoid. 

One way to make progress is to use the parameter set 
determined from the susceptibility data and to exam- 
ine whether other properties can be understood with the 
same parameter set as well. Here we choose to study the 
specific heat. In the upper panel of Fig. [8] the phonon 
frequencies are assumed to be known a priori; so the free 
phonon contribution added is the one for this known fre- 
quency u> = loq. Clearly, none of the static spin models 
describes the susceptibility data and the specific data sat- 
isfactorily. From the knowledge of both quantities com- 
pelling evidence can be deduced that a dynamic spin- 
phonon coupling must be present. 




FIG. 8: Comparison of the specific heat for the parameter sets 
determined from the susceptibility in Fig. |H] For the upper 
plot the phonon frequency is assumed to be known before- 
hand. For the lower plot no a priori knowledge is assumed for 
the phonon frequency; so it is also fitted. The line styles are 
the same as those used in Fig. [7] that means the exchange 
couplings used are those from the corresponding curves in 
Fig.0 



But the situation can be less advantageous in practice. 
Let us assume that we do not possess knowledge about 
the frequency of the phonon to which the spin system 
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may or may not be coupled. Then this frequency could 
be seen as an additional fit parameter. This view point 
was adopted in the lower plot in Fig. |H1 The parame- 
ter set (a = 0.1; 5 = 0) of the dashed-dotted curve can 
be clearly discarded. The parameter set (a = 0; S = 0, 
the Heisenberg model HM) of the solid curve can be dis- 
carded because it does not describe the susceptibility, cf. 
Fig. □ The parameter set (a = 0.12; S = 0.05) of the 
dashed curve seems to fit the data down to T an 0.7 J. 
For lower values, however, the agreement is poor. So 
also this data set must be discarded. Note that we refer 
only to the temperature regime where the extrapolated 
series yield reliable results. 

So we have shown for the above example that it is not 
possible to describe the data of a dynamic spin-phonon 
system in the intermediate regime, where all energy scales 
are of similar magnitude, in the framework of static spin 
models plus (decoupled) phonons. But it is not sufficient 
to consider only one quantity at (relatively) high tem- 
peratures. In order to obtain unambiguous evidence of 
the presence of a dynamic spin-phonon coupling one has 
to dispose of either data down to low temperatures (not 
considered here) or to consider at least two independent 
thermodynamic quantities. Otherwise, one may easily 
be misled by a good agreement in one quantity alone to 
conclude that a simpler static model is sufficient for a 
microscopic description. 

IV. SUMMARY AND OUTLOOK 

The problem of a one-dimensional spin system coupled 
to phononic degrees of freedom is investigated. A clus- 
ter expansion is applied to obtain a series expansion in 
the magnetic coupling J. Results are computed at finite 
temperatures for the free energy, the specific heat, and 
magnetic susceptibility. No truncation in the phononic 
subspace is necessary since the expansion is performed in 
J about the limit J = 0, not in the inverse temperature. 
This is the first realization of a cluster expansion at finite 
temperatures for an extended spin-phonon problem. The 
implementation of the expansion in a computer program 
yields high orders in the expansion parameter. 

The comparison of the results obtained from the ex- 
trapolated series to data from quantum Monte-Carlo sim- 



ulations shows a very good agreement. This supports the 
reliability of the approach. Hence our results can serve 
as input for quantitative data analysis since the features 
of the magnetic susceptibility and the specific heat at 
moderate and at high temperatures T > 0.15 — 0.25 J are 
described reliably. 

A possible route to improve the extrapolations is to 
understand the limit J — > oo in the Hamilton operator 
better. At first glance, this limit looks simple since 
it corresponds to the adiabatic situation with lu/J — ► 0. 
But this limit is not straightforward since the nature of 
the excitations is unclear at present. Note that in this 
limit static displacements can be made at no energetic 
cost. This suggests that the excitations are domain walls. 
Whether these objects are static (because to /J — > 0) or 
dynamic (because the magnetic subsystems retains its 
fluctuations) is unclear at present and constitutes an in- 
teresting theoretical issue in itself. 

Finally, we analyzed the data of a spin-phonon system 
in the intermediate coupling regime, where all energies 
are of similar magnitude, in great detail. As expected, we 
could show that the effects of the dynamic spin-phonon 
coupling cannot be imitated by a static spin model and 
decoupled phonons. But it is necessary to study low tem- 
peratures or at least two quantities like the susceptibility 
and the specific heat carefully. Otherwise, one may be 
easily misled by a good agreement in one quantity alone 
to conclude that a static microscopic model is sufficient. 
We think that this conclusion is helpful for future exper- 
imental analyses of low-dimensional spin systems. 
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